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I INTRODUCTION 


In  an  earlier  paper  (Ref.  [1]),  equations  of  motion  for  the  transient 
response,  to  a shock  wave,  of  a submerged  shell  with  internally  attached 


1 

• li 


structures  were  obtained  by  a substructuring  procedure.  Several  quantities 
in  these  equations,  relating  to  the  attachments,  were  only  broadly  defined, 
and  determined  in  detail  only  in  an  ad  hoc  manner  for  specific  appendages. 

In  this  report,  matrix  expressions  for  these  quantities  are  defined  which 
permit  their  computation  for  arbitrary  internal  structures  approximated  by 
finite  elements. 

The  present  report  addresses  the  same  basic  problem  as  Ref.  [1],  but 
does  not  consider  the  details  of  the  fluid  loading  on  the  shell,  which  may  be 
found  in  Ref.  [1].  The  effects  of  external  loads  applied  to  the  internal 
structures,  which  were  not  considered  in  Ref.  [1],  are  included  here. 

The  substructuring  procedure  described  herein  is  intended  for  use  in 
transient  response  problems.  Unlike  many  other  substructuring  schemes,  which 
are  used  primarily  to  solve  steady-state  problems  (see,  e.g..  Ref.  [2]),  the 
present  method  does  not  involve  the  modes  and  natural  frequencies  of  the  combined 
system  (shell  and  appendages).  In  addition,  a system  stiffness  matrix  is  not 
required,  since  the  present  technique  employs  the  interaction  forces  and 
moments  at  the  junctions  of  the  shell  and  the  attached  structures. 
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II  GENERAL  KINEMATICS 


Consider  a shell  S of  arbitrary  geometry  (see  Fig.  1).  Let  n,  2 and  3 

denote  mutually  perpendicular  directions,  respectively  normal  and  tangent  to 

the  middle  surface  of  the  shell.  If  d , d„  and  d-  denote  corresponding 

n z j 

displacement  components  of  this  surface,  a 3 by  1 displacement  matrix  d^(s,  t) 
may  be  defined,  in  which  s denotes  an  arbitrary  point  on  the  middle  surface. 


t is  time,  and  underlining  indicates  a matrix;  d is  assumed  positive  inward. 

n 


Corresponding  to  this  displacement  matrix,  the  3 by  1 jth  shell  mode 

has  components  (ji^,  (j)^  and  (})j,  with  which  are  associated  a natural  frequency 


and  a generalized  mass 


md) . (b  dA 


in  which  m denotes  the  mass  per  unit  of  surface  area  A of  the  middle  surface 
of  the  shell  and  the  superscript  T denotes  a transpose.  The  shell  modes, 
which  satisfy  the  orthogonality  condition 


m dA  = 0 if  k ^ j 


may  be  used  as  expansion  functions  for  the  shell  displacements,  i.e., 

v_ 


d(s,  t)  = I q„.  (t)  ^ (s) 
k=l  ^ 


in  which  the  q (t)  are  generalized  coordinates  for  the  shell  and  V is 

O K o 


the  number  of  modes  of  the  shell  considered. 


(1) 


(2) 


(3) 


Assume  that  an  elastic  structure  O,  entirely  within  the  shell  S,  is  attached 
to  S at  points  Sj(I  = 1,  2,...,  N^) , as  shown  in  Fig.  2.  Let  O be  approximated 
by  finite  elements  with  interface  (attached)  nodal  points  I,  coinciding  with 


I 

ei 


Sj,  and  interior  (unattached)  nodal  points  a(a  = 1,  2,...,Ny),  corresponding 


to  material  points  of  o not  In  contact  with  S.  The  finite  element  representation 
of  O is  assumed  to  include  rotational  degrees  of  freedom. 

At  each  Interface  point  I,  certain  components  of  displacement  and/or 
rotation  of  o must  match  corresponding  quantities  of  S,  depending  on  the 
type  of  support.  Thus,  at  an  interface  point  I,  define  the  vector 

( 

X 

-ol 


where  the  upper  partition,  an  N^.  by  1 matrix,  contains  those  physical  degrees 

of  freedom  of  O which  must  match  corresponding  physical  degrees  of  freedom  of 

/\ 

S,  and  the  lower  partition,  an  by  1 matrix,  contains  the  remaining  physical 

*) 

degrees  of  freedom  of  a active  at  interface  point  I. 


Since  the  interior  points  of  o are  not  connected  to  S,  define  the  vector 
of  physical  degrees  of  freedom  of  0 at  an  interior  point  a as 


X 

—00, 


(5) 


It  may  be  observed  in  Eqs.  (4)  and  (5),  and  in  what  follows,  that  the  symbol  - 
above  a quantity  indicates  a quantity  associated  with  physical  degrees  of 
freedom  of  o constrained  to  move  with  S,  while  the  symbol  ^ above  a quantity 
indicates  a quantity  associated  with  physical  degrees  of  freedom  of  o not 


At  a simple  support,  e.g.,  the  displacements  of  O and  S must  be  equal, 
but  the  rotations  of  a and  S are  not  equal.  Thus,  assuming  two  active 
physical  degrees  of  freedom,  the  displacement  of  o is  located  in  the 
upper  partition  of  Eq.  (4),  and  the  rotation  of  o is  located  in  the 
lower  partition  of  Eq.  (4). 


- A - 


constrained  to  move  with  S. 


If  the  appropriate  definitions  of  Eqs.  (4)  and  (5)  are  applied  to  all 
nodal  points  of  0,  the  physical  degrees  of  freedom  of  0 may  be  arranged  in 
the  form 


X 

-o 


X 

—a 


(6) 


in  which  the  upper  partition,  an  N by  1 matrix,  contains  the  terms  from 
all  interface  points  [Eq.  (4)]  and  the  lower  partition,  an  N by  1 matrix, 
contains  the  terms  from  all  interface  points  [Eq.  (4)]  and  the 
terms  from  all  interior  points  [Eq.  (5)]. 

In  the  sequel,  the  elements  in  the  upper  partition  of  Eq.  (6)  will  be 
referred  to  as  "constrained  physical  degrees  of  freedom  of  o",  while  the 
elements  in  the  lower  partition  will  be  referred  to  as  "unconstrained  physical 
degrees  of  freedom  of  o".  The  matrix  ^ has  size  N by  1,  where 


N = N + N 


(7) 


is  the  total  number  of  active  physical  degrees  of  freedom  of  a. 


From  Eq.  (3),  the  shell  displacements  at  an  interface  point  1 may  be 
written  as 

v„ 


d(s^,  t)  = d^  = I qg^(t)  i^(Sj)  i 
k=l 


(8) 


where  is  a 3 by  Vg  matrix  with  elements  and  g[g  is  a Vg  by  1 matrix 

containing  the  generalized  coordinates  of  the  shell.  Similarly,  the  rotations 
of  the  shell  cross  section  at  an  interface  point  1 may  be  expressed  as 


k»l 


(9) 
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in  which  denotes  the  3 by  1 matrix  of  rotations  about  the  n,  2 and  3 

axes  in  the  mode  ^ and  is  a 3 by  v matrix  with  elements  4)^(s  ). 

Continuity  of  appropriate  components  of  displacement  and  rotation  of 
shell  and  substructure  at  an  Interface  point  I requires  that 


(10) 


where  D^,  an  by  6 matrix,  is  a local  rotational  coordinate  transformation 
which  transforms  quantities  expressed  in  the  S-coordlnate  system  into  the 

o-coordinate  system.  Each  row  of  contains  at  most  three  non-zero  terms. 

Use  of  Eqs.  (8)  and  (9)  permits  Eq.  (10)  to  be  written  in  the  form 

(11) 

(12) 

Application  of  Eq,  (11)  to  all  interface  points,  with  reference  to  Eq.  (6) 
for  proper  ordering,  permits  the  continuity  between  shell  and  substructure 
to  be  written  as 

^ = (13) 


where  is  an  N^.  by  matrix  given  by 


in  which  C is  an  N by  matrix  having  the  form 


The  coordinate  transformation  array  in  Eq.  (10)  is  used  for  convenience 
in  presentation.  Actually,  only  a 3 by  3 coordinate  transformation  array 
is  required  at  each  Interface  point,  the  rows  of  this  array  applying  to 
both  displacements  and  rotations. 


(lA) 


Equation  (13)  expresses  the  constrained  physical  degrees  of  freedom  of  a in 

*) 

terms  of  the  free-free  modes  of  S. 

As  suggested  in  Ref.  [3],  the  response  of  a may  be  taken  as  the  super- 
position of  its  static  response  to  the  actual  motion  of  its  supports  and 
the  dynamic  response  with  respect  to  fixed  supports,  i.e.**\ 

As  /\ 

In  Eq.  (15),  is  the  N by  fixed  modal  matrix,  where  denotes  the 
number  of  fixed-base  modes,  and  is  the  by  1 matrix  of  corresponding 

_ As 

generalized  coordinates.  The  N by  N matrix  ^ is  obtained  by  considering 
the  static  response  of  O to  an  arbitrary  support  motion  (support  "settlement") 
when  the  nodal  loads  corresponding  to  the  unconstrained  physical  degrees  of 
freedom  vanish.  In  such  a case,  Eq.  (15)  becomes 

/V 

^^^static  ^ ^^^static 

*) 

Equation  (13)  replaces  Eq . (5)  of  Ref.  [1]: 

^(Sj,  t)  = d(Sj,  t) 

**) 

Equation  (15)  replaces  Eq.  (A)  of  Ref.  [1]: 

^(s,  t)  = ^ ^1^®^  ■’■gCs.  t) 
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As  shown  in  Ref.  [2],  the  matrix  g (sometimes  termed  the  "constraint 
modes")  may  be  obtained  by  considering  the  linear  static  "force-displacement" 
relation  for  O: 


= R 

o static  — 


(17) 


in  which  K is  the  N by  N unconstrained  stiffness  matrix  for  the  substructure 
and  R is  the  corresponding  vector  of  nodal  forces  and  moments.  For  the 
determination  of  Eq.  (17)  may  be  partitioned  to  correspond  to  Eq.  (6), 
yielding 


” j 

-21  i 

-12 

X 

-o 

" I " 

i^22 

1 

0 

— 

L _J 

Static 

— - 

(18) 


since  the  only  nodal  loads  are  the  reactions.  It  then  follows  from  the  lower 
partition  of  Eq.  (18)  that 


(X  ) 

-a  static 


-k"^  k 
-22  -21 


(X  ) 

“O  static 


(19) 


where  the  superscript  -1  indicates  an  inverse.  A comparison  of  Eqs.  (19) 
and  (16)  shows  that 


g ~ "^22  -21 


(20) 


Thus,  the  elements  of  g are  the  static  displacements  and  rotations  (including 
rigid-body  motion),  corresponding  to  the  unconstrained  physical  degrees  of 
freedom  of  a,  which  result  from  successive  unit  displacements  and  rotations, 
corresponding  to  the  constrained  physical  degrees  of  freedom  of  a,  at 
attachment  points  I when  all  other  constrained  physical  degrees  of  freedom 
of  a are  set  equal  to  zero. 


The  use  of  Eq.  (13)  permits  Eq.  (15)  to  be  written  as 


Xj  = io  So  + £ C gg 


(21) 
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As  suggested  In  Ref.  [2],  Eqs.  (13)  and  (21)  may  be  combined  Into  the  single 


matrix  equation 


X = T q 
-o  - ^ 


• • c 

in  which  is  given  by  Eq . (6),  the  transformation  matrix  T has  the  form 


C I 0 


I ^ 


and  the  array  of  generalized  coordinates  is  given  by 


In  Eq.  (22),  5^  is  an  N by  1 matrix,  £ is  a V by  1 matrix,  with 


and  T is  an  N by  V matrix.  Equation  (22)  may  be  written  in  subscript 


notation  as 


"ak  - J,  ’j  ^ « 

J = 1 
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III  EQUATIONS  OF  MOTION  OF  SHELL 


The  kinetic  and  potential  energies  of  the  combination  of  shell  S and 
internal  structure  0 may  be  written  as 


, S N 

T = i I y.qo-+T  I M.  x^. 

^ ^2  S3  2 j Oj 


j=l 

v„ 


N N 


1 o 9 « , 111 

" = 2 " 2 ‘^jk  ^ok 


(27) 


(28) 


in  which  the  are  given  by  Eq.  (26),  the  are  the  elements  of  the 

partitioned  stiffness  matrix  of  Eqs.  (17)  and  (18),  and  the  are  the 

elements  of  the  partitioned  mass  matrix  of  a from  an  assumed  lumped  mass 
*) 

formulation  . The  mass  matrix  must  be  partitioned  to  correspond  to  Eq.  (6), 
i.e. , 


M 


ei 

1 

0 

0 1 

in^ 

— 1 

-22 

(29) 


where  the  upper  left  partition,  an  N by  N matrix,  and  the  lower  right 

/V  ^ 

partition,  an  N by  N matrix,  are  diagonal. 


The  equations  of  motion  of  the  shell  may  be  derived  from  Lagrange's 
equations : 

dt^3qc5j^  3^  ° ^Sj  " l,2,...,Vg) 


(30) 


*) 


If  Eq.  (26)  is  substituted  into  Eqs.  (27)  and  (28),  one  obtains  the 
equations  which  replace  Eqs.  (6)  and  (7)  of  Ref.  [1]: 


10 


in  which  the  are  the  generalized  forces  corresponding  to  the  total  fluid 
pressure  p acting  on  the  surface  of  the  shell  (see  Fig.  1)  and  the  Q , are 
the  generalized  forces  from  the  action  of  0 on  S,  exclusive  of  the  contribution 
derivable  from  the  kinetic  and  potential  energies.  As  shown  in  the  Appendix, 

/s 

the  Q . are  the  elements  of  the  matrix 
^ J 

^ T ^ 

Ss  = C 1 R 

where  R,  an  N by  1 matrix,  is  the  vector  of  known  applied  nodal  loads  acting 
on  a (see  Fig.  2) . 


It  follows  from  Eqs.  (26)  and  (24)  that 


9X  , 
ok 

9"ak 

(j  = 1,  2,. 

T,  ■ 

k = 1,  2,. 

Vg; 


which  may  be  seen  to  reference  terms  in  the  upper  and  lower  left  partitions 
of  Eq.  (23).  Thus,  the  use  of  Eqs.  (27),  (28),  and  (32)  leads  to 


dt^9qgj^  " ^Sj  \ ^ak 

N N 

9V  2 V c V 

5 = U.to.  q^.  + ) T,  . y K,  . X . 

3 .1  kj  ki  01 


The  substitution  of  Eqs.  (31),  (33),  and  (34)  into  Eq . (30),  with  reference 
to  Eqs.  (18),  (23),  and  (29),  results  in 


''j  Vi  ^ “j‘"i  '’Sj  = ^Tj  - 2 Vg) 


where  the  P are  the  elements  of  the  vector 
J 


P = P + P 
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whose  addends  are  defined  as 


T^T 


tji  V !Si2 

(37) 

+ h2i  + !<22  i - SI 

(38) 

The  bracketed  term  in  Eq.  (38)  may  be  recognized  as  the  physical 
equations  of  motion  of  0.  Hence,  the  contribution  of  Eq.  (38)  vanishes, 
and  the  equations  of  motion  of  S become 


Uq  +yu)  Qn, . Q . 

j^Sj  j ^Sj  ^Tj  ^oj 


(39) 


in  which  the  generalized  forces  , corresponding  to  the  forces  and  moments 


exerted  by  0 on  S,  are 


Oaj-- j,  j,  'S.i  V 


k=l 


i=l 


(40) 


which  result  from  the  dynamic  reactions 


R 


Ok  ■ >‘ak  * \l  *ol> 


exerted  by  0 on  S 


*) 


(41) 


Using  Eqs.  (22)  to  (24),  the  equations  of  motion  of  S may  now  be  written 
in  the  form 


*) 


Equations  (40)  and  (41)  are  generalizations  of  specific  examples 
given  by  Eqs.  (17),  (18);  (22),  (23);  (Al) , (A2) ; and  (A12),  (A13), 
(A16)  of  Ref.  [1]. 
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N 

k=i  'k“jk  ■ .1^  ■ "j-j  ^SJ 


(42) 


N 

«TJ  - '^IJ^  <J  - '•  2 V 


where  6.,  is  the  Kronecker  delta  and  the  F . , the  elastic  forces  and  moments 
Jk  i 

exerted  by  a on  S and  vice-versa,  are  the  elements  of  the  N by  1 matrix 


l - (Js-u  + iiij  £)  £ % + hi  i % («) 

The  effect  of  inertia  forces  appears  in  the  double  summation  in  Eq . (42) 

It  may  be  observed  that  the  equations  of  motion  of  the  shell  have  stiffness 
coupling  only  between  shell  and  substructure. 

In  the  equations  of  motion  of  the  shell,  the  contributions  of  the  Inertia 
forces  and  the  elastic  reactions  of  a on  S may  be  given  a simple  physical 
interpretation  by  referring  to  Fig.  3.  Figure  3a  shows  an  idealized  finite 
element  model  of  a substructure  O attached  to  a shell  S,  and  Figs.  3b  and  3c 
show  the  manner  in  which  Eqs.  (42)  and  (43)  treat  the  connection  of  o to  S. 

The  effect  of  inertia  forces  results  from  the  rigidly  attached  point  mass  of 
Fig.  3b,  while  the  effect  of  the  elastic  reactions  of  O on  S results  from  the 
"spring  force"  of  Fig.  3c. 


The  equations  of  motion  of  the  substructure  O may  be  obtained  from 
Lagrange's  equations: 

(:  = 1.  2....,  v^) 

Oj  OJ 

in  which  the  are  the  generalized  forces  corresponding  to  the  known 
applied  nodal  loads  acting  on  a.  As  shown  in  the  Appendix,  the  Q^. 
are  the  elements  of  the  matrix 


It  follows  from  Eqs.  (26)  and  (24)  that 

^ ^^gk  ^ c (j  = 1,  2,.. ; 

^^gj  k = l,  2,...,N) 


which  may  be  seen  to  reference  terms  in  the  upper  and  lower  right  partitions 
of  Eq.  (23).  Since  the  upper  right  partition  of  Eq.  (23)  is  a null  matrix, 
the  use  of  Eqs.  (27),  (28),  and  (46)  leads  to 


d 9T 

^ak+N 


dv  _ r X V 

^‘igj  ^ k=l  i=l  ""gi 

The  substitution  of  Eqs.  (45),  (47),  and  (48)  into  Eq.  (44)  results  in 


^Okj^\+N  ^gk+N  ‘Sc+N,!  ^gi  “ ° 


which  may  be  written  in  matrix  form  as 


^ foi  S22I  Jg_  i!^22^ 


(50) 
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If  Eqs.  (22)  to  (24)  are  substituted  Into  Eq.  (50),  and  the  matrix  products 
are  carried  out,  there  results 

£ ^22  ^ 5 ^22  ^ + 5(^21  + ii22  ^ % 


+ 4 1522  £ C ig  = ^ R 


In  view  of  Eq.  (20),  which  '^efines  the  constraint  modes,  the  third 
term  of  Eq.  (51)  vanishes.  Moreover,  as  discussed  in  Ref.  [2],  the  fixed- 
base  modes  of  a are  the  solutions  of  the  eigenvalue  problem 


i22%>J  ■ “oi  =22 


where  0)^^  is  the  natural  frequency  of  the  jth  mode,  and  hence  are  orthogonal. 
Thus,  Eq.  (51)  may  be  simplified  to  yield  the  equations  of  motion  of  a : 


where  and  are  diagonal  matrices  whose  main  diagonal  terms  are, 
respectively,  the  generalized  masses  and  generalized  stiffnesses  of  a,  and 


r’  -1  IT 
G = ^ 5,22  £ 


is  a by  N matrix  whose  columns  are  the  expansion  coefficients  of  a modal 

series  representation  of  the  constraint  modes.  In  subscript  notation,  Eq.  (53) 
*) 


becomes 


* “»J  "oj  ■ Vj  \ - =jk  ' <”) 


It  may  be  observed  that  the  equations  of  motion  of  the  substructure  a have 
Inertia  coupling  only  between  shell  and  substructure. 


Equation  (55)  may  be  recognized  as  a generalization  of  Eq.  (A8)  of  Ref.  [1]. 
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V CONCLUDING  REMARKS 


The  equations  of  motion  for  the  transient  response  of  a shell  S with 
internal  structure  o are  given  by  Eqs.  (A2) , (43),  (53),  and  (54).  It  is 
important  to  note  that  no  total  system  (S  + o)  modal  stiffness  matrix  is 
required  in  the  present  formulation,  as  is  required  when  eigenvalues  and 
eigenvectors  for  the  entire  system  are  sought  (see,  e.g,.  Ref.  [2]).  It 
should  also  be  noted  that,  in  view  of  the  stiffness  coupling  between  S and  a 
in  Eqs.  (42)  and  (43),  the  modal  mass  matrix  to  be  Inverted  in  the  solution 
of  a transient  response  problem  has  size  v by  V , determined  solely  by  the 
number  of  modes  of  the  shell  employed  in  the  solution.  Formulations  with 
inertia  coupling  between  shell  and  substructure  require  the  size  of  the  modal 
mass  matrix  to  be  inverted  to  increase  as  the  total  number  of  modes  v increases. 


The  physical  stiffness  and  mass  matrices  of  a,  K and  M of  Eqs.  (18)  and 
(29),  respectively,  may  be  obtained  by  modifying  a finite  element  computer 
code  such  as  SAPIV  (Ref.  [4]).  The  constraint  modes  may  be  determined  by 
means  of  Eq.  (20),  or  by  utilizing  a solution  option  available  with  SAPIV.  Hence, 
it  is  possible  to  evaluate  all  of  the  matrices  required  in  Eqs.  (42)  and  (43) 
by  using  only  those  rows  of  K and  only  those  main  diagonal  terms  of  M which 
correspond  to  the  constrained  physical  degrees  of  freedom  of  a.  Once  the 
expansion  coefficients  of  the  constraint  modes  [Eq.  (54)]  and  the  modal 
coefficients  for  the  interaction  forces  and  moments  [Eq.  (43)]  have  been 
evaluated,  the  constraint  modes  and  the  matrices  K and  M may  be  eliminated  from 
the  final  form  of  the  equations  of  motion  of  S and  a.  This  elimination 
and  the  absence  of  a system  modal  stiffness  matrix  will  reduce  the 
computer  core  (central  memory)  required  to  solve  a given  problem. 


kikuiiattH 
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The  present  formulation  will  also  reduce  the  amount  of  computer 
time  (central  processor  time)  required  to  set  up  the  matrices  in  the  equations 
of  motion,  as  compared  with  that  required  by  a method  employing  a system  modal 
stiffness  matrix.  By  way  of  example,  if  the  formulation  of  Ref.  [2]  is  used 
for  a transient  response  problem,  the  matrices 


,j,CT 

K T^ 

(■^6) 

CT 

c 

I 

M T 

(57) 

must  be  evaluated.  The  computer  time  required  to  evaluate  Eq.  (56)  may  be 
excessive,  especially  if  the  problem  under  consideration  requires  many  modes 
V and/or  many  physical  degrees  of  freedom  of  a,  and  if  K is  not  well  banded. 

It  should  be  pointed  out  that  the  rearrangement  of  physical  degrees  of  freedom 
of  a to  produce  the  partitioning  of  Eq . (18)  usually  results  in  a poorly  banded 
matrix.  For  a lumped  mass  formulation,  the  computer  time  required  to  evaluate 
Eq.  (57)  is  not  significant  compared  with  that  needed  to  evaluate  Eq . (56). 

Thus,  a reduction  in  the  amount  of  computer  time  required  to  set  up  the 
equations  of  motion  is  obtained  by  the  elimination  of  a system  modal 
stiffness  matrix.  A reduction  in  solution  time  is  also  quite  probable. 

The  free-free  modes  of  the  shell  S may  be  determined  by  means  of  a computer 
code  such  as  B0S0R4  (Ref.  [5]),  while  SAPIV  may,  of  course,  be  used  to  obtain 
the  fixed-base  modes  of  the  substructure  a . Tf  a consistent  mass  formulation 
is  desired  for  a,  the  equations  of  Sections  III  and  IV  may  be  modified  in  a 
straightforward  manner.  However,  in  such  a case,  SAPIV,  which  is  based  on 
a lumped  mass  formulation,  must  be  replaced  by  a computer  rode  employing  a 
consistent  mass  formulation.  Although  the  equations  derived  in  Sections  III 
and  IV  deal  with  one  piece  of  internal  structure,  the  present  method  is  easily 
extended  to  Include  any  number  of  internal  appendages. 
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APPENDIX  - GENERALIZED  FORCES  DUE  TO  APPLIED  LOADS  ON  SUBSTRUCTURE 

The  virtual  work  due  to  the  known  applied  nodal  loads  on  substructure  o 
may  be  written  as 

5W  = (6]^)  R 

in  which  6X  is  the  N by  1 matrix  of  virtual  displacements  and  R is  the 
-a  — 

/s 

N by  1 matrix  of  applied  loads.  In  Ea . (Al) , both  the  virtual  displacements 
and  the  applied  loads  correspond  to  the  unconsti alned  physical  degrees  of 
freedom  of  0 . From  Eq . (21),  it  follows  that 


so  that  Eq . (Al)  may  be  expressed  as 


rn  Arp  A 'T'  T*  ^ 'P  ^ 

6W  = 63^ (^  R)  + 6£g(c'  s,  R) 


The  virtual  work  of  Eq.  (Al)  may  be  written  in  terms  of  generalized 
forces  as 


6W  = 6^  5 


where  is  given  by  Eq . (24)  and 


s = 


-^s 


So 


is  a V by  1 matrix  partitioned  to  correspond  to  Eq.  (24) 
Eqs.  (24)  and  (A5)  in  Eq.  (A4)  results  in 

6W  = 6qJ 


The  use  of 


(Al) 


(A2) 


(A3) 


(A4) 


(A5) 


(Ah) 
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Since  the  variations  of  the  generalized  coordinates  are  arbitrary,  a 
comparison  of  Eqs.  (A3)  and  (A6)  yields 

''  T ^ 

° - I - (A7) 

% " K - (A8) 

The  components  of  generalized  force  associated  with  the  generalized  coordinates 
of  the  shell  S are  given  by  Eq . (A7) , and  those  associated  with  the 
generalized  coordinates  of  the  substructure  o are  given  by  Eq.  (A8) . 
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